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ABSTRACT 

The natural history of the biogenic elements (H,C,0,N) from their first association 
within cold molecular clouds to their delivery to the Earth during the late bombardment 
of the inner solar system, is intimately linked to water ice. The earliest organic 
compounds are formed in cold interstellar molecular clouds as a result of UV and 
thermal processing of sub-pm ice grains which contain trapped carbon and nitrogen 
molecules. Structural changes in the water ice host underlie and fundamentally control 
important macroscopic phenomena such as the outgassing of volatiles, the rates of 
chemical reactions, and processing and retention of organic compounds. Prebiotic 
organic material was in all likelihood delivered the early Earth in a pristine state as a 
consequence of its sequestration within a protective water ice host. 


INTRODUCTION 

Most scientists believe that terrestrial life as we know it originated on the Earth. 
However, observations made over the past 30 years suggest that many of the organic 
compounds antecedent to life could not have formed on the pre-biotic Earth. The 
early Earth is now thought to have had a more oxidizing atmosphere of C0 2 , rather 
than the CH 4 atmosphere postulated earlier. The primitive oceans, rather than being a 
“primordial soup” were in all likelihood much more dilute, incapable of sustaining 
chemical reactions of sufficient variety or quantity for pre-biotic organic syntheses. 
Scientists now believe that much of the early prebiotic organic chemistry in fact took 
place in the extraterrestrial environment - within ices in the cold molecular cloud 
which collapsed to form the early solar nebula, and in icy and rocky planetesimals 
(the source material of comets and asteroids/meteorites, respectively) of the outer 
solar nebula. 

Three sets of observations provide evidence that organic compounds exist in great 
quantity in the extraterrestrial environment. First, analyses of the Murchison and 
Allende carbonaceous chondrite meteorites, which fell to Earth in 1969, reveal 
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abundant organic compounds. A list of these compounds is shown in Table 1. The 
parent bodies of these meteorites, which are among the most primitive and ancient 
objects in the solar system, must have either accumulated from materials containing 
these organic compounds, or the compounds themselves must have been synthesized 
in the parent bodies early in their history. Second, remote telescopic observations in 
the infrared and radio wavelengths show that simple organic compounds exist in 
abundance within cold molecular clouds (Table 1 ). Third, spectral analyses of chemical 
species present in the gaseous outflow of comets reveal a variety of simple organic 
species (Table 1 ). These species must have been frozen as simple ices in the cometary 
nucleus which volatilized upon warming. 

THE ORIGIN OF THE BIOGENIC ELEMENTS AND THEIR FIRST 
ASSOCIATION WITHIN INTERSTELLAR ICE 

OVERVIEW 

Mixed ices containing mostly water are primary components of cold interstellar 
molecular clouds which are the birthplace of stars and protoplanetary nebulae (see 
Figure 1). This interstellar ice represents the earliest and most primitive association 
of the biogenic elements. Frozen on sub-micron-sized silicate grains, water ice serves 
as a matrix for the low-temperature ( 1 5-50 K) trapping of a variety of small carbon- 
bearing molecules. Exposure to UV radiation forms highly reactive free radicals 
which upon warming recombine to create a refractory organic residue. Laboratory 
simulations of the growth, radiation processing and thermal annealing of astrophysical 
ice analogs show that a variety of organic molecules can be synthesized, including 
complex hydrocarbons, alcohols, amines, nitrides, and esters or ketones as well as 
organic acids (Table 2). Remote infrared observations demonstrate that these organic 
compounds can exist in the interior of interstellar molecular clouds for long periods 
of time at ultra-low temperatures ( 1 0 - 20 K), protected from high particle or radiation 
fluxes by the bulk of the cloud itself. 

Most low melting point solids were volatilized during gravitational collapse within 
the cold molecular cloud which formed our solar nebula. However, in the outer solar 
nebula, original icy condensates from the parent cold molecular cloud may have 
persisted or may have been volatilized and recondensed to form new ices representative 
of the early outer solar nebula. These condensates coalesced to form icy planetesimals 
(the parent bodies of comets) which now reside in the Oort cloud and the Kuiper belt. 
Early in solar system history, some of these icy planetesimals bombarded the volatile- 
poor rocky planets of the inner solar system, and thus provided the Earth with its rich 
complement of gases, volatiles, water, and biogenic elements and compounds. 

Cometary bodies (Figure 2) are the most pristine and primitive objects in the 
solar system, possibly containing relict biogenic compounds and frozen volatiles from 
the parent molecular cloud. Thus, if one could deduce the physical and chemical 
state of components and phases within a cometary nucleus as it exists and can be 
observed today, insight could be gained into the nature of the volatiles and biogenic 
compounds introduced to the Earth ~4 billion years before the present. 


Table la 


Concentrations and Molecular Characteristics of Soluble Organic Compounds 

in Meteorites [73] 


Class 

Concentration 

(PPM) 

# ComDounds 
identified 

Chain Lenath 

Amino Acids 

60 

74 

c 2 -c 7 

Aliphatic hydrocarbons 

>35 

140 

c r C 23 

Aromatic hydrocarbons 

15-28 

87 

c 6 -c 20 

Carboxylic acids 

>300 

20 

C 2 _C 12 

Dicarboxylic acids 

>30 

17 

c 2 -c 9 

Hydrocarboxylic acids 

15 

7 

c 2 -c 6 

Purines and Pyrimidines 

1.3 

5 

NA 

Basic N-heterocycles 

7 

32 

NA 

Amines 

8 

10 

c,-c 4 

Amides 

55-70 

>2 

NA 

Alcohols 

11 

8 

c r c 4 

Aldehydes & Ketones 

27 

9 

c r c 5 

Total 

>560 

411 



Table 1b 


Organics Detected in Astrophysical Ices [73-76] 



Cometary* 

Interstellar Ice 

Molecule 

Abundance 

Abundance 

h 2 o 

100 

100 

CO (in polar ice) 

5-10 

1-10 

CO (in non-polar ice) 

~ 

10-40 

CH 3 OH 

1-5 

<4-10 

C0 2 

3 

1-10 

XCN 

.1 

1-10 

H 2 

-1 

-1 

ch 4 

0.2-4. 5 

-1 

HCO 


-1 

h 2 co 

0-5 

1-4 

nh 3 

0 . 1-2 

1-4 

N 2 

- 0.02 

10-40 

0 2 

10-40 

10-40 

0 3 

few 

few 

OCS or CO 3 

few 

few 


Molar concentrations relative to water ice 


Fig. 1. Gas Pillars in Ml 6 
the Eagle Nebula. The 
dark pillar-like structures 
are columns of cool 
interstellar hydrogen gas 
and dust that are 
incubators for new stars. 
Inside cold molelcular 
clouds such as these, 
water-rich ices are 
condensed as thin rinds 
onto sub-gm dust grains. 
These interstellar ices 
represent the most 
primitive association of the 
biogenic elements, and 
are the sites of the first 
organic syntheses. Credit: 
Jeff Hester and Paul 
Scowan (ASU) and NASA. 



Much of our present knowledge of the composition of interstellar ices is due to 
detailed laboratory studies of the growth and thermal and UV processing of mixed 
ices. These studies consist primarily of fingerprint-style matching of astronomical 
IR spectra with laboratory IR spectra from astrophysical ice analogs. The physics of 
cometary behavior has been studied by IR and gas release studies of analog ices as 
well as by direct observation of cometary behavior in the IR and radio wavelengths. 
This work has allowed scientists to infer the composition of various icy bodies and/or 
to explain on a phenomenological basis the outgassing events which are associated 
with coma formation. 

However, a knowledge of the petrology (structure, phase relationships, 
microstructural associations and assemblages) of cometary and astrophysical ice 
analogs is fundamentally important because such physical properties underlie and 
control chemical reaction rates, volatile retention and release, vaporization behavior, 
thermal conductivity, bulk density, infrared spectral characteristics and so on. For 
solid phases it is really the structure of the icy materials which we would like to 
determine. Once structures are known, one can predict, correlate or explain in a 
fundamental way remote observations of real astronomical objects. In the research 
described herein, the detailed structure of crystalline and amorphous ices is determined 
directly by electron diffraction and related to other laboratory and astronomical 
observations. 


ICE IN THE INTERSTELLAR MEDIUM 

Cold molecular clouds, the birthplace of stars and planetary nebulae, are dynamic 
environments in which small molecules exist both in the gas phase and as solid phases 
(principally water-rich ices) condensed onto small silicate and carbon grains. Gas 
phase molecules are observable through milimeter and sub-milimeter radio wavelength 
observations [1-3]. Using interferometry techniques, images with a resolution of 
1 00- 1 000 AU (one astronomical unit is the distance from the Sun to the Earth) can be 
obtained which allow the investigation of temperature and gas composition within 
some star forming regions. Solid phase molecules, usually in the form of water-rich 
ice mantles on grains, are observable through infrared (IR) absorption measurements 
against a background illumination source, such as a star embedded within a cold 
molecular cloud [4-6]. 

Gas phase abundances can be modelled using gas phase ion-molecule reaction 
schemes. However, these reactions are found to be insufficient to explain the observed 
data of some species, indicating that small molecules are continually replenished 
from the solid state reservoir (i.e., from volatile-rich icy grains). The occurrence and 
abundance of gas phase molecules, in combination with short accretion timescales in 
cold molecular clouds, imply that a fast (-105 year) gas phase replenishment occurs 
[7,8], This replenishment involves outgassing of ices and desorption of molecules, 
both determined by the structure of the solid host which is usually H 2 O. A great deal 
has been inferred about the bulk composition and chemistry of ice in molecular clouds 
from laboratory studies of astrophysical ice analogs [5,9-19]. These laboratory studies 
are continuing and expanding in response to the availability of recent IR astronomical 
observations of clouds as well as data from the Infrared Space Observatory. 

The emphasis has traditionally been on IR spectral data which can be used for 
“fingerprint” matching of observed absorption bands in cold molecular clouds. 
However, IR spectroscopy does not have the capability to deduce structural or phase 
information from materials in the solid state, although typically such information is 



inferred indirectly to interpret the data. It is necessary therefore to utilize a structure- 
sensitive technique (i.e., X-ray, neutron or electron diffraction) in laboratory studies 
of astrophysical ice analogs in order to understand or deduce the mechanisms 
underlying the observed phenomena. 

THE INTRODUCTION OF WATER ICE, ORGANIC MOLECULES AND 
VOLATILES INTO THE EARLY SOLAR NEBULA AND TO THE EARTH 

Most low melting point solids such as water ice were volatilized during cloud 
collapse and the formation of the solar nebula. However, in the outer solar nebula, 
original icy condensates from the parent cold molecular cloud either persisted 
unchanged, or volatilized and recondensed to form new ices representative of the 
early outer solar nebula. These condensates coalesced to form icy planetesimals 
which became a part of the Oort cloud and the Kuiper belt. The Oort cloud is a 
spherical region of space extending from 1000 - 30,000 AU, beginning far beyond 
Pluto and continuing to halfway to the nearest star. The region is believed to contain 
millions of icy objects 1-50 km diameter, in highly elliptical orbit about the sun. 
These objects were formed near the present orbit of Jupiter, and cast out of the inner 
solar sytem by Jupiter’s gravity [20], The Kuiper belt is a ring-shaped region just 
outside the orbit of Pluto also believed to contain millions of icy objects. Kuiper belt 
objects, unlike Oort cloud objects, are thought to have remained at the same distance 
from the sun since their formation [20]. 

The inner rocky planets of the solar system, Mercury, Venus, Earth and Mars, are 
believed to have been devoid of volatiles at the time of formation of the Oort cloud 


and the Kuiper belt - due to the high temperatures experienced at this distance from 
the proto-Sun. However, early in solar system history, a large number of icy 
planetesimals bombarded the volatile-poor rocky planets of the inner solar system, 
and in all likelihood provided the Earth with its rich complement of gases, volatiles, 
water, and biogenic elements and compounds. The evidence for this late bombardment, 
a period of intense cratering early in solar system history, had been largely 
circumstantial until lunar rocks were dated, and the cratering histories of the Moon, 
Mars, and moons of other planets were assembled. The Earth, having a plate tectonics 
cycle as well as a cycle of weathering, has retained almost no record of this early 
bombardment. 


THE IMPORTANCE OF STRUCTURAL DATA FOR EXPLAINING 
WATER ICE BEHAVIOR 

LABORATORY STUDY OF ASTROPHYSICAL ICES 

The study of volatile release from amorphous and crystalline water ice has direct 
application to cometary and planetary outgassing phenomena as well as to studies of 
the thermal properties of astrophysical ices. Several investigations report anomalous 
gas release or vaporization phenomena which are as yet only poorly understood. For 
example, [2 1 , 22] have described 8 different temperature regimes of gas release in 
water ice as it is warmed from 16° K to melting, yet water undergoes only 2 phase 
transitions in this temperature range. Two (possibly three) of the temperature regimes 
of gas release which occur below the amorphous to crystalline transformation suggest 
that changes in morphology, porosity, or short range local order (“annealing”) are 
occurring even in the amorphous state. Sandford and coauthors [23] report different 
physical trapping mechanisms for different H 20 :C 0 ratios in their ices based on a 
detailed analysis of the fine structure of infrared spectra. Kouchi [24] studied the 
vapor pressure of amoqrhous and crystalline ices, and found an apparent two order of 
magnitude increase in vapor pressure for amorphous ice relative to crystalline ice at 
the same temperature. These data were used to explain the gas production rate of 
new comets vs. old. Kouchi and Kuroda [25, 26] studied the gas release and crystalline 
structure of low temperature mixed ices. Kochan and coauthors[27] report results 
from the KOSI experiments, in which large amounts of cometary ice analogs are 
grown and processed by heat and irradiation to release gases. Hsuing and Roessler 
[28] in their analyses of these simulated cometary nuclei from the KOSI experiments, 
report phase separation and migration of frozen volatiles within the ice. However, 
the original ices for the KOSI experiments were formed by spraying water and dust 
into liquid nitrogen. Under these conditions, the resulting ice is in its high temperature 
hexagonal form and the results are not directly applicable. Other studies [29-31] 
have been performed which address the presence or absence of ice clathrates in solar 
system objects and in laboratory analogs. The presence or absence of clathrates bears 
strongly on the physico-chemical nature of the ices and their structural and gas release 
properties. However the only direct proof of the presence of an ice clathrate is by 
analysis of structural data contained in diffraction patterns. 



STRUCTURAL STUDIES OF ASTROPHYSICAL ICE ANALOGS 


Nearly all of the laboratory studies of astrophysical ice analogs carried out to 
date record secondary or tertiary manifestations of structural or morphological change 
rather than observing such changes directly. While this approach is useful in mimicking 
the bulk physical and chemical phenomena taking place in cometary and other 
extraterrestrial ices, it does not reveal the mechanisms which control such phenomena. 
This phenomenological interpretation of data is probably the cause of somewhat 
contradictory explanations invoked to account for differences in behavior in similar 
temperature regimes. Furthermore, even in the case of [25, 26], structural data obtained 
by reflection electron diffraction are integrated over large volumes of ice relative to 
the size of individual crystallites. Thus, these structural studies do not reveal crystalline 
imperfections and cannot be used to characterize microstructural relationships, grain 
size, crystallite orientation and shape, grain boundary effects and inter- and 
intracrystalline defects. It is this microstructural, micromorphological and 
microchemical heterogeneity which must be characterized if the mechanisms 
underlying the observed phenomena are to be understood. 

In order to elucidate the mechanisms responsible for the phenomenological 
observations described above, direct determinations of the structure and composition 
of astrophysical ice analogs prepared under realistic conditions are required. These 
analyses should be obtained at or below the level of resolution of individual crystallites 
within the ices (<5-100 nm). It is highly desirable to image individual domains of 
short-range order within amorphous ices (if they are present), to observe nucleation 
and growth phenomena from the vapor state onto grain surfaces, to follow structural 
changes during temperature-induced annealing, and to view the ices at high resolution 
during entrapped gas release phenomena. Likewise, it is desirable to study, at high 
resolution, the structural state, morphology, and microchemistry of mixed ices below, 
at, and above all of the temperatures at which subtle IR spectral features (attributed to 
different local ice “architecture;” see [23]) appear or change. 


ELECTRON MICROSCOPY OF WATER ICE 

LOW-PRESSURE FORMS OF WATER ICE 

Water ice at low pressure can exist in three metastable amorphous forms, one 
metastable crystalline polymorph, and one stable crystalline polymorph. Amorphous 
water ice was first described early in this century. Diffraction studies were first 
performed by Burton and Oliver [32], Low-density amorphous water ice (Iaj) occurs 
when water vapor is deposited slowly onto an inert surface at 77 K. Iai has a density 
of 0.94 [33-35] and has been extensively studied because of its similarity to liquid 
water [36-39], and because it can be used as an embedding medium for biological 
tissues observed in electron microscopy [40-42], High-density amorphous water ice 


(Iah) was first reported by [34], X-ray diffraction patterns of water vapor deposited 
onto a substrate cooled to 10 K showed the presence of oxygen atoms at distances 
between those of the first and second nearest neighbor at the interstitial site of the Iai 
structure, which results in a higher value for the calculated density of oxygen atoms 
(i.e., 1.1 instead of 0.94). For this reason, Narten and coauthors [34] introduced the 
names “high density amorphous ice”, (Iah) and “low density amorphous ice” (Iaj). 
Subsequent experiments failed to reproduce the high-density form and from this it 
was concluded that the high density amorphous form was an artifact of epitaxial 
growth on an oriented single crystal [37, 43-44], Later it was shown that high density 
amorphous ice can form during vapor deposition on amorphous carbon films [41, 
45], Another form of high density amorphous water ice can be made by pressure 
melting, and is also referred to as high-density amorphous ice (hda or HDA). When 
cooled to 77 K, crystalline hexagonal ice, Ih, undergoes a transition into hda at 1-2 
GPa pressure [46, 47] . The resulting ice has a density of 1 . 1 7 and can persist metastably 
for some time at low pressure and 77 K. When this pressure- induced high density 
amorphous ice is warmed the ice transforms irreversibly into low density amorphous 
ice in the range 1 05- 1 23 K [46, 48-5 1 ]. The radial distribution function of hda [49,52] 
differs significantly from that of the Iah of [34] and it has been argued that pressure- 
induced hda has more hexagonal order than vapor-deposited Iah [53]. 

When amorphous water ice is warmed to a high enough temperature, the ice 
crystallizes into the cubic crystalline polymorph (Ic, space group Fd3m) before, at 
higher temperatures, recrystallization to the stable hexagonal crystalline polymorph 
(Ih, space group P 63 /mmc) occurs. Prior to cubic crystallization, amorphous water 
ice undergoes a weak glass transition that manifests itself as an endothermic step in 
differential thermal analysis (DTA) scans [54-57], During this transition, the 
amorphous solid transforms into a viscous liquid. During the crystallization process, 
a significant amount of amorphous water ice persists as a third amorphous form which 
has a slightly different structure than Iaj [33,45]. This form of amorphous water ice 
is called the “restrained amorphous form” (Ia r ) because it apparently inhibits 
crystallization into the cubic crystalline polymorph. 

The cubic crystalline polymorph (Ic) can result from warming amorphous ice 
through the amorphous - cubic ice transition temperature ( 1 20- 1 40 K, time-dependent), 
from direct deposition of water vapor onto a substrate in the temperature range 1 30- 
1 50 K, or from direct deposition of water vapor onto a substrate at temperatures < 1 20 
K, due to heat of fusion, if deposition rates are rapid enough. In fact, crystallization 
to Ic (as a consequence of heat of fusion) during vapor-deposition at T< 1 20 K occurs 
in all but the slowest deposition experiments, with the result that making bulk quantities 
of low-density amorphous water ice in the laboratory is not possible. 

The hexagonal crystalline polymorph results form the cooling of liquid water 
below the freezing point, from the vapor-depositon of water above -160 K, and from 
the warming of cubic crystalline ice above -230 K. The hexagonal crystalline 
polymorph is the most common form of water ice on the Earth, manifesting itself in 
the familiar form of snowflakes. 



ELECTRON MICROSCOPY OF ASTROPHYSICAL ICE ANALOGS 


A brief description of the modifications required to deposit, image and process 
astrophysical ice analogs in an Electron Microscope is given in [58]. The instrument 
is shown in Figure 3. Modifications were made to the microscope to accomplish two 
purposes: First, vacuum system modifications were made to reduce the partial pressure 
of water vapor and hydrocarbons by several orders of magnitude so that vacuum- 
born sample contamination is substantially reduced. Second, a low electron dose 
system was fabricated and installed with the intent of minimizing damage to electron 
beam sensitive samples. Images are collected digitally using a cooled, slow scan 
CCD camera and processed using a commercial image processing program. A gas 
inlet tube was installed so that gases of selected compositions could be leaked into 
the sample chamber and vapor-deposited on a cooled sample substrate. AQuadrupole 
Mass Spectrometer was interfaced directly into the objective lens polepiece and is 
used to monitor gas deposition conditions, and to record pressure changes for up to 
12 masses during the course of an experiment. 

Samples are introduced into the electron microscope in two ways: (1). Ices 
collected or grown elsewhere can be inserted into the microscope in the frozen state 
using a cryotransfer device. The liquid nitrogen stage used for cryotransfer operates 
in the range 100-373 K. It has a cryoshutter which allows the sample and specimen 
rod to be taken out of the microscope while at low temperature. Frozen samples can 
be loaded and unloaded from the stage inside a cryotransfer station and transported at 
liquid nitrogen temperatures in special dewars. (2). Ices can be grown in situ by 
leaking pre-mixed gases directly onto a cold sample substrate inside the microscope. 
The specially fabricated liquid helium cooled stage is equipped with a shutter to 
protect the specimen during extended periods of annealing. A heater control allows 
temperatures between 1 5 and 373 K. 

The substrates for ice deposition are ~5 nm thick carbon films supported by 400 
mesh Cu electron microscope grids. A typical in situ deposition and annealing 
experiment proceeds as follows: A sample bulb filled with water vapor or other gas 
mixtures at reduced pressure ( 1 -20 Torr) is placed on the microscope, and the sample 
insertion lines are evacuated. The cryostage is then put in the microscope and cooled 
to the desired starting temperature. While the sample substrate is being viewed at 
low magnification, the gas is leaked into the chamber surrounding the substrate surface. 
A needle valve is used to control the flow of gas into the sample chamber and the 
resulting vapor deposit of ice is monitored by inspection. Ice can be deposited at the 
rate of nanometers per hour to microns per hour. The sample substrate can be tilted 
directly toward the gas inlet (for ballistic deposition), or oblique to the inlet so that 
molecules must diffuse indirectly to the substrate. The latter method of deposition 
disallows the formation of water clusters due to supersonic flow, as described by 
[59], 

Figure 4a shows an electron diffraction pattern from the ice which indicates that 
it is amorphous at the level of electron diffraction (some short-range order is certainly 
present, but no trace of crystallinity is evident in the pattern). The sample is warmed 
up slowly and diffraction patterns are taken at regular intervals of temperature. Figure 



Fig. 3. Electron Microscope, modified for high vacuum and low electron-dose imagery of 
ices. (1). Gas bulb used to leak pre-mixed gases onto the sample substrate inside the 
microscope. (2). Gas inlet valve, which can be used to adjust the rate of deposition of ices 
from nanometers per hour to microns per hour. (3). Cryogenic sample holder, capable of 
temperatures from 15-363 K. (4). Quadrupole mass spectrometer used to record mass 
spectra from gases released during thermal annealing. (5). Ion pump which, in concert 
with a LN 2 anticontamination device, maintain a low-10 7 torr vacuum at the sample. (6). 
Digital camera with image intensifierfor recording low-dose images of ice. (7). TV monitor 
for real-time viewing of images and diffraction patterns. (8). Computer for recording and 
analyzing digital images and diffraction patterns in real time. 

4b shows the first indication of the amorphous to cubic transition (Ia r — >Ic). The 
arrow in the figure marks the first indication of an I c 220 maximum. This transition 
occurs at 143° ± 2° K, which is (within error) the transition temperature reported in 
the literature [40], However, as we will show later, the temperature of onset of cubic 
crystallization is a sensitive function of annealing time and the rate of warming of the 
ice, and may vary by as much as 20-30 degrees, depending on temperature history. 
Figures 4c and 4d show further increases in crystallinity. 

In addition to viewing electron diffraction patterns and recording them 
electronically, it is also possible to record moderately high resolution images at low 
electron dose in order to view the morphology and defect structure of individual 
crystals and aggregates at nanometer resolution. In crystalline samples, it is possible 
to view the size, shape, and degree of perfection of individual crystallites by using 
bright and dark field imaging. Figure 5a shows a selected area electron diffraction 



pattern for hexagonal ice obtained from water vapor deposited at 160° K. Figure 5b 
shows a cubic diffraction pattern, obtained from water vapor deposited at 140° K. By 
tilting the electron beam so that different parts of the diffraction pattern are allowed 
to pass through the objective aperture, bright and dark field images are formed. Figure 
5c for example, is a bright field image of cubic ice, obtained by directing the central 
beam through the aperture (position labeled as “a” in Figure 5b). Figure 5d shows a 
dark field image from the same area, obtained by placing the aperture in the position 
labeled “b” in Figure 5b. By careful dark field tilting experiments, it is possible to 
identify the types of defects in individual ice crystallites, as well as to determine their 
orientation, size and shape. 


REAL-TIME ELECTRON DIFFRACTION AND COMPUTER REDUCTION 
OF PATTERNS 

The phase changes and structural annealing which take place in amorphous and 
crystalline water-rich ices are critically dependent on time and temperature. In gas 
release and IR studies of ice, the rate of data collection is typically not a limiting 
factor in an experiment (that is, spectra can be recorded rapidly enough to identify 
changes as they occur). However, for conventional analysis using X-ray diffractometry, 
a single high resolution pattern may require hours of intergration time. In the classic 


Fig. 4. Electron 
diffraction patterns of 
low-temperature water 
ice. (a). Diffraction 
pattern of a 10pm 
diameter region of 
amorpnous ice which 
was vapor deposited at 
100 K and warmed to 
133 K. Note the two 
diffuse maxima 

(arrowed) but no sharp 
maxima indicative of 
crystalline ice. (b). 
diffraction pattern of 
amorphous ice recorded 
at 1 43 K at the onset of 
the amorphous to 
crystalline cubic ice 
transition. The arrow 
shows a sharp band of 
intensity at the position 
of the 220 maximum of 
cubic ice. (c). Diffraction pattern recorded at 1 53 K showing the increase in crystallization 
of the cubic phase, (d). diffraction pattern recorded at 163 K showing further transformation 
to cubic ice. It is noteworthy that the amorphous ice never fully transforms to the crystalline 
state until complete recrystallization into the hexagonal polymorph. 




work of Narten and coauthors [34] for example, ice samples were deposited over a 
period of 2 weeks, and structural interpretations were based on only 4 diffraction 
patterns. In our experiments the temperature is increased at the rate of 1 K per minute 
(analogous to typical IR experiments), and it would be desirable to be able to collect 
patterns at the rate of at least one per Kelvin, over temperature ranges of 200 K. By 
conventional film methods, an electron diffraction pattern can be exposed in as little 
as 2 seconds, but the microscope only carries a small number of sheets of film per 
cassette. This quantity is insufficient particularly when one wants to bracket exposures 
and collect standardized patterns before and after the experiment is run. 

Hardware and software improvements to the imaging system were made so that 
diffraction patterns could be collected on-line in real time. As a result, hundreds of 
patterns can be collected at the rate of one every few seconds, and stored electronically. 
Computer programs reduce the two-dimensional diffraction patterns (figure 6a) to 
one-dimensional radial intensity functions (figure 6b). These intensity vs. radial 
distance functions are calibrated using an internal crystalline gold standard which is 
sputtered onto a small portion of the support substrate, and an amorphous carbon 
background is subtracted to represent the 5 nm substrate. Then, an atomic scattering 


background is subtracted to yield a diffracted intensity distribution in reciprocal space 
(figure 6c). 


LOW-TEMPERATURE STRUCTURAL TRANSFORMATIONS IN WATER ICE 


Figure 7 shows background-subtracted radial intensity distributions of ice 
diffraction patterns during warming from 1 5 to about 200 K. High density amorphous 
ice (Iah) forms when water is vapor deposited below about 65 K. This high density, 
low-temperature form of amorphous ice transforms irreversibly to low density 
amorphous ice (Iai) above about 65 K. Low density amorphous ice gradually 
transforms to cubic ice at 140-150 K. If one subtracts a pure crystalline cubic ice 
pattern from a crystallized amorphous ice pattern, the residual intensity resides in a 
third amorphous phase (Ia r ), which persists to higher temperatures (>190K) along 
with the crystalline ice. It is notable that no more than about 30% of the amorphous 
ice ever crystallizes. This has important implications for the volatilization of water 
and trapped gases from comets. 

These patterns are striking, and reveal much about the dynamics of water ice. 
However, to interpret real-time data recorded using other techniques (IR, gas release, 
etc. ), one needs to interpret the data more quantitatively. Computer software was 






written which fits Gaussian curves beneath the two principal maxima in the amorphous 
water ice pattern. In this way, one can study small changes in the position, height and 
breadth of the peaks as a function of temperature or processing history. Two parameters 
of the amorphous diffraction maxima were found which reflect repeatable and 
substantive changes in the structure of the amorphous ice during warming. Figure 
8a-b shows these parameters for a series of about 1 50 patterns recorded at 1 per K per 
minute during an experiment. Each data point in the figure represents a measurement 
of a particular feature of a particular diffraction maximum at a specified temperature. 
From these data one can quantify the transformation from one amorphous form to 
another and relate these changes to other real time changes. For example, the position 
of the first amorphous peak appears to faithfully reflect the transition of high density 
amorphous ice to low density amorphous ice (Iah — > Iai), which occurs between about 
35 and 65 K. Changes in the intensity of the first amorphous peak reflect the 
transformation from the low density amorphous form to the restrained amorphous 
form (Iai — > Ia r ) between 130 and 150 K. 

The Iah -> Iai transformation between 35 and 65 K can be correlated with and 
indeed, underlies and explains other data, for example changes in the character of CO 
trapping sites [23], changes in the shape of the 3.07 pm water band as a function of 
temperature [60], or the temperatures at which radicals recombine in UV photolyzed 
ices to form organic compounds [6 1 -62], These and other conclusions are more fully 
explained in [45, 63]. These papers conclude that high density amorphous ice (Iah) is 
the most common form of ice in the universe, and that the Iah — > I a l transformation 



Fig. 7. Electron diffraction 
radial intensity distributions 
(diffractograms) of the 
different forms of vapor- 
deposited water ice from a 
region about 10 pm in 
diameter. The forms 
encountered are high- 
density amorphous ice (Iah), 
low-density amorphous ice 
(Iai), restrained amorphous 
ice (Iar) (while transforming 
from Iai), and a composite of 
restrained amorphous and 
cubic crystalline ice (Ic). The 
latter can be separated into 
two components by 
subtraction of a standard 
cubic pattern as indicated, 
which allows one to follow 
changes in the amorphous 
component after the onset of 
cubic crystallization (after 
[45]). 



controls the recombination of biogenic elements to form organic compounds within 
interstellar ices. This research provides a physical basis for subtle changes in the 
3.07 pm water band in astrophysical ice analogs and provides ground truth for the 
determination of the temperature and to some extent the thermal history of 
astrophysical ices. Small changes in the IR band structure of molecules trapped in 
the ice (e.g., CO), can be used to illuminate the structure state and thermal history of 
the parent ice. The Iah -» Iai transformation has important implications for ices 
within cold molecular clouds as well as for ices within comets. 


THE STRUCTURE OE AMORPHOUS WATER ICE; EXPERIMENTAL DATA 
AND MOLECULAR DYNAMICS SIMULATIONS 


Fourier transformation of the diffraction patterns of amorphous water ice yields 
the radial distribution function (rdf) of distances between oxygen atom pairs [64, 65]. 
Figure 9 shows rdfs from our electron diffraction data along with X-ray diffraction 
data from pressure amorphosed ice [49], and X-ray diffraction data from vapor- 
deposited ice [34], The high density form is characterized in all cases by a splitting 
of the second nearest neighbor peak such that one component lies between the first 
and second nearest neighbor peaks of the low-density amorphous form. Molecular 
dynamics calculations were made of the formation of amorphous water ice (details of 




Fig. 8a. Temperature depencence of the position of the first and second diffraction maxima 
in the temperature range in which amorphous water ice transforms from the high-density 
form to the low-density form. The transition is well defined, but sluggish. The ice was 
vapor-deposited at 1 5 K and held for 37 minutes before warm-up at a rate of 2 K per minute 
(after [45]). 
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Fig. 8b. Temperature depence during warmup of the intensity (in arbitraty units) of, from 
top to bottom, the second and first diffraction peaks of the amorphous component and the 
cubic 220 peak, which marks the onset of crystallization into Ic. The onset of the transition 
la, -> la r is well defined by an increase in intensity of the first diffraction peak. The heating 
rate was 1 to 2 K per minute, and the onset of the cubic crystallization occurred 182 min. 
after depositon at 15 K (after [45]). 
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this are shown in [63, 66]) and these results were compared with experimental data. 
The similarity in the rdf results suggests that the modeling correctly represents the 
ice deposition process. Once one is reasonably certain that the computer model of 
ice is a good stand-in for the real ice formed in the electron microscope experiment, 
the computer generated structure can be sampled to determine the average coordination 
number of the oxygens, the strength and number of hydrogen bonds per atom, the 
size of pore spaces in the structure, etc (e.g., see Table 3). The data appear to support 
the model of [50] in which the high density amorphous form is a collapsed low- 
density lattice, characterized by the presence of oxygen-oxygen distances between 3 
and 4 A. 

Jenniskens and Blake [45] propose the existence of two fundamentally different 
types of cometary ice based on thermal history. Two families of comets are known: 
Kuiper belt objects, which remain during their lifetime at temperatures below 30 to 
50 K, and Oort cloud objects, which were expelled from their formational zones in 
the Uranus-Neptune region where the ice was heated to some 60 to 100 K. After 
being expelled, Oort cloud comets cool down to temperatures T<20 K but retain the 
structural signature of their higher temperature origin. Kuiper belt objects should be 



made of high density amorphous ice, whereas Oort cloud comets should be made of 
low density amorphous ice. The difference in ice structure may manifest itself in the 
outgassing characteristics of new comets. 


HIGHER TEMPERATURE STRUCTURAL TRANSFORMATIONS IN 
AMORPHOUS WATER ICE AND ICE CRYSTALLIZATION 

A second amorphous-to-amorphous transformation is seen in vapor-deposited 
amorphous ice (Iaj) as it is heated to a temperature just below to the crystallization 
temperature (Figure 8b). This transformation results in a new form of amorphous 
water ice called the restrained amorphous form (Ia r ). The restrained amorphous form, 
as its name implies, resists crystallization, and indeed, persists metastably with cubic 
crystalline ice throughout its range. The onset of the Iai -> Ia r transition coincides 
with the glass transition (T g ) in water ice, a structural relaxation which marks the 
change from a solid to a viscous liquid [67, 68]. It is this transition, and not the 
crystallization of amorphous ice per se, which underlies the release of trapped gases 
and initiates many other processes (e.g., clathrate formation in impure ices; [69]) 
which occur in the temperature range 120-150 K. During the Iai —> Ia r transition, 
pore spaces are annealed out of the structure and residual dangling -OH bonds are 
eliminated. This results in the release of “excess volatiles” held in pore spaces, while 
the restrained amorphous ice is still capable of retaining at least 10 % volatile species 
until it volatilizes or is recrystallizes into cubic or hexagonal ice. 

Structural data from time-temperature-transformation experiments indicate that 
no more than 30% of amorphous water ice ever crystallizes to the cubic polymorph 
below about 200 K [33, 45]. This resistance of Ia r to crystallization is examined in 
[70] , where the first direct evidence is presented that amorphous water ice above its 
glass transition behaves as a “strong” liquid (Angell, 1993). Nucleation rates of 
cubic crystalline water ice within a Ia r precursor are unlike those measured for liquid 
water clusters [71], suggesting that restrained amorphous ice is thermodynamically 
distinct from liquid water. The temperature dependence in the model that describes 
the liquid water clusters is too steep below T = 150 K to fit the nucleation rates for Ia r . 
The behavior of Ia r as a strong liquid manifests itself in the persistence of the 
amorphous form to high temperatures and its coexistence with cubic crystalline ice 
over long time periods. This anomalous behavior suggests that within cometary ices, 
the bulk of the ice volatilizes directly from the amorphous state rather than first 
crystallizing to cubic or hexagonal ice. 

[72] pointed out that the “strong liquid” character of amorphous water ice above 
the glass transition can be modified by the incorporation of impurities. It is likely 
that Ia r will have a differing viscosity-temperature relationships based on the nature 
and quantity of impurities in the ice. It is interesting to note that, if a more “fragile” 
liquid character is imparted to the amorphous ice due to impurities, the steep 
dependence of viscosity on temperature will result in the persistence of amorphous 
ice to higher temperatures or for longer periods of time. 

Lastly, we note that there must be some structural cause for the “strong liquid” 



character of Ia r . What prevents Ia r from fully transforming into cubic ice? The 
answer may be that Ia r contains a recognizable short-range hexagonal stacking order 
which restrains the materials from crystallizing [45]. To form cubic crystals, the ice 
has to fully restructure all layers with ABAB (hexagonal) type stacking into ABCABC 
(cubic) type stacking. Domains with hexagonal short-range order probably resist 
restructuring into cubic short-range order and persist metastably until wholesale 
restructuring transforms all ice into the thermodynamically stable hexagonal 
polymorph at higher temperatures. Electron micrographs of amorphous water ice 
lend some insight. Figure 1 0 shows Ia r coexisting with cubic ice after various annealing 
times. The inset diffraction pattern shows some hexagonal ice (Ih) maxima despite 
the fact that the ice is in the cubic temperature domain. Many of the crystallites in the 
image show planar defects, often separated by a few tens of Angstroms, which we 
interpret to be stacking faults. These stacking faults are probably mistakes in ABCABC 




Fig. 10. Bright field image of water ice, vapor-deposited in the amorphous state and 
annealed to 132 K for 24 hours. Note individual crystals of ice, some showing planar 
defects with spacings less than 5 nm (arrowed). The defects are attributed to local 
stacking disorder in which the normal ABCABC... stacking sequence of the cubic 
polymorph is interrupted by hexagonal ABAB... planes. 
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ordering, which would result in small hexagonal domains in the ice. 
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